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Abstract. A direct numerical solution of the steady-state Boltzmann equation in a cylindrical 
geometry is reported. Finite-size effects are investigated in large semiconducting nanowires using 
the relaxation-time approximation. A nanowire is modelled as a combination of an interior with 
local transport parameters identical to those in the bulk, and a finite surface region across whose 
width the carrier density decays radially to zero. The roughness of the surface is incorporated by 
using lower relaxation-times there than in the interior. 

An argument supported by our numerical results challenges a commonly used zero-width 
parametrization of the surface layer [Q. In the non-degenerate limit, appropriate for moderately 
doped semiconductors, a finite surface width model does produce a positive longitudinal magneto- 
conductance, in agreement with existing theory [Q. However, the effect is seen to be quite small (a 
few per cent) for realistic values of the wire parameters even at the highest practical magnetic fields. 
Physical insights emerging from the results are discussed. 



1. Introduction 

> ■ 
• i— i . 

^ 1 The effect of a finite system size on the conductivity of a material is a subject of considerable 
cd ' physical interest, which has recently been lent added relevance and importance by rapid develop- 
ments in nanowire synthesis and assembly El HI 13 03 El El E3 EI] , electrical characterization and 
transport measurement methods [HI C2J CHI EH CHI CHI CZl CHI CHI E01 EH |22] - The development 
of nanowires currently represents an important part of materials and applied physics research. Hi- 
erarchial self-assembly techniques ^3] envisaged in semiconductor nanowires make them promising 
central elements of future integrated electronics. 



Reports of basic functional two- and three-terminal semiconductor nanowire devices including 
junctions, bipolar transistors and field-effect transistors are now widely found in the literature. The 
nano-scale transport properties of several important semiconductors including Si, EH E21 E| , 
GaAs 1213, GaN[IH] and Ge [231 EHE5] ) , and semi-metals EZl EHl EHl EDI EU E21 ESJ have been 
investigated in detail. (Bi has attracted great attention because of its unique combination of inter- 
esting properties and its potential for thermoelectric applications [2B1 I2E1 EH ED E21 ESI EH EHj) 
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Bicrystalline nanowires crossed nanowire structures [12], functional networks [12112111121 includ- 
ing ultra-high-density lattices [T7j, heterostructures [HE] and superlattice devices [2U| are part of the 
rapidly growing body of novel nanowire configurations under development. 

A central aspect of theoretical enquiry must be the extent to which the conductivity of a 
nanowire differs from that of the bulk material. At first glance, it seems reasonable to suppose that 
the conductivity is smaller in a nanowire because of the addition of surface scattering, assuming 
that the band structure does not change drastically. However, experimental results have yielded 
conflicting indications on this point, which is presently not well-understood. Particularly intriguing 
are reports of mobility values higher than their bulk counterparts observed in silicon nanowires |18j : 
while in other cases the mobility has been deduced to be orders of magnitude lower. 

Such challenging theoretical questions, brought into immediate relevance by the extensive data 
on electrical transport in nanowires compiled in the last few years, make a thorough quantitative 
investigation valuable at this point. The widespread pursuit of experiments pertinent to the surface 
effect on conductivity motivates a generic numerical description of the finite-size effect allowing 
both freedom and simplicity in the incorporation of nanowire characteristics. 

For large- and moderate- sized nanowires operated at room temperature, semi-classical kinetic 
effects are expected to be as important as quantum mechanical effects like the modification of band 
structure. Thus a semi-classical study is both necessary and desirable in the common regime where 
the diameter of the nanowires is much larger than the thermal wavelength. While a quantum me- 
chanical approach is indispensable to investigate conductance in very narrow nanowires (radii of a 
few nm) |37| E?5] , it is neither viable nor suitable for larger nanowires with radii of the order of 100 
nm. For the latter the use of the semi-classical Boltzmann equation is the appropriate method to 
approach a thorough quantitative study of finite-size effects. 

These effects have previously been addressed theoretically only with simplistic assumptions 
regarding the surface. The analytical results most widely quoted are due to Chambers P, who 
used kinetic-theoretical ideas to calculate the modification of the effective mean free path due to 
the presence of a surface and then used this value to find the conductivity. Essentially the same 
results were obtained using the Boltzmann equation with suitable boundary conditions in 

We present here the results obtained from a general numerical scheme we have developed to 
solve the Boltzmann equation in a cylindrical geometry. Such a direct numerical solution offers the 
requisite freedom in incorporating nanowire characteristics and conditions such as the equilibrium 
electron density profile, or the presence of defects and impurities. To retain conceptual simplicity, 
we employ the relaxation-time approximation as a first step towards the systematic modelling and 
prediction of nanowire conductivity with a view to relating model parameters to experimental data, 
and possibly casting a light on surface characterization [HI EU]- I n this context, comprehensive 
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experimental studies of Bi nanowires are of direct relevance ["""""I I3*H 13*21 13*3*1 . The introduction of 
diameter-controlled synthesis of nanowires [El E] provides yet another fruitful context for our study. 

We emphasize that surface scattering is just one of many scattering mechanisms that contribute 
to the resistivity of a nanowire; other mechanisms, especially acoustic phonon scattering, can be 
produce more dissipation in many circumstances. In addition, for studies of specific nanowires with 
specific surface defects and impurities, the relaxation time approximation employed here would ide- 
ally be replaced by a more exact approach such as Monte Carlo or quantum mechanical simulation 
as mentioned above. Our intention, however, is to contribute generic intuition about the behavior 
of surface scattering in large wires, even when the specific surface scattering centers are unknown 
and even when other resistive effects may be primary. For this reason, we use a generic relaxation 
time approach and map out its predictions for varied choices of wire characteristics. Despite its 
approximate character, this approach confers useful insight and is substantially more refined than 
the often-invoked "specularity coefficient" model |TJ[?""""|- 

Inter alia, the results yield insights into the general problem of transport with a nominal re- 
laxation time that varies with spatial coordinates, which is non-trivial because of the fact that the 
diffusion of carriers connects different spatial regions, making their properties inter-dependent, and 
thereby introduces a connectivity to the physical situation. This may be of direct relevance to 
transport in layered media such as magnetic multi-layers [3U [321 HBH 133] • 

The numerical problem resulting from a finite-difference representation of the Boltzmann equa- 
tion with a simple grid is solved using the conjugate- gradient method, which offers significant com- 
putational efficiency. A comparison with analytical results in limiting cases confirms the reliability 
of the scheme. 

The paper is organized as follows: In section [21 we describe our computational framework, 
and the form of the Boltzmann equation adapted to the problem at hand. Our surface model is 
described in section El We present our results in the simple case of zero-magnetic-field in section 
EJ showing how they conform to physical expectation, thus validating the numerical scheme used. 
The important limit of zero-surface width is considered in section [""] Section [fj] is devoted to the 
longitudinal magneto-conductance arising due to the surface effect and is followed by a summary 
(section [TJ) . 

2. Boltzmann Equation and Problem Specification 

In our approach, we consider a cylindrical conductor which has a finite surface width, so that there 
is no abrupt change at the boundary to be dealt with through boundary conditions. In particular, 
this includes treating the unperturbed (i.e. when the external fields are zero) distribution function 
fo as a function of both space and momenta. In order to use the Boltzmann equation with such 
a distribution function f we introduce an effective internal electric field E c (r) in addition to any 
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external field, to account for this spatial variation of fo. In a confined cylindrical system where fo 
decays from its interior magnitude to zero over a finite surface width, the force due to this internal 
field is exactly analogous to the constraining force that keeps a particle within an enclosure with 
an abrupt boundary, and is therefore physically expected. With this additional internal field, the 
Boltzmann equation jU] for the distribution /(p, r) of non- interacting carriers in the relaxation 
time approximation with a spatially varying relaxation time is 

V/ • — + Vp/ • q(E + E c + — x B) = -Lz*! (1) 
m m r(r) 

Here E and B are the external electric field and magnetic field, and q and m are the charge and 
mass of each carrier. It is convenient to introduce the deviation 0(p, r) = /(p, r) — /o(p, r) due to 
the presence of the external fields; the equation for <ft is 

4> = -r(r)[V0 ■ - + qV P <P ■ (E + E c + — x B) + qV P f • (E + — x B)]. (2) 

m m m 

We consider a long cylinder with no azimuthal or axial inhomogeneity and confine attention to the 
case where the external fields are uniform and parallel to the axis of the wire, defining a natural 
axis z of reference: E = Ez; B = B z (see Fig. QJ) 

Thus the physical problem at hand requires only one spatial degree of freedom r = ^/x 1 + y' 2 , 
although we have to treat all the momentum components. It is convenient to work with the local 
momentum components p r = p • f and pg = p • 9 rather than the canonical momenta in a cylin- 
drical system, because the canonical momentum conjugate to 9 is an angular momentum. Also, to 
preserve form (2) of the Boltzmann equation, /, / and are defined as coordinate densities in the 
space of (p r ,pg,p z ,r) so that fdp r dpgdp z dr is the number of carriers in a volume element of this 
space. The missing volume element factor 2nr is absorbed into /, changing the internal field E c , 
as we shall see. We solve equation (2) by discretizing on a real space grid of about 20 4 grid points. 
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The resulting matrix equation is solved by the method of conjugate gradients. 

The drift velocity, which is the population average of the component of the carrier velocity 
parallel to the external electric field, is 

/ v ) = = 1 f d 3 pdrf(p,r)p z = 1 J d 3 p dr<p(p,r)p z 
m m fd 3 pdrf(p,r) m f d 3 p drf (p,r) 

The effective conductivity a is q(v z )N/lA where N is the total number of carriers in the wire, I is 
its length and A its cross-sectional area. 

An explicit expression for E c can be directly deduced from the spatial variation of /o for an 
arbitrary form of the latter, since by definition of the unperturbed distribution, 

V/ • - + Vp/o • <?E C = (4) 

m 

For simplicity, we assume that the spatial dependence of the distribution function is completely 
separable from the momentum dependence so that fo = F(p)£(r) and the internal field Ec points 
only in the direction of r : Ec = Ec{r)f. The computational results described below pertain to 
the case where the momentum-distribution is Maxwellian: 

Inserting fo = F(p)£(r) and (5) in (4) we see that the simplest consistent form of the internal field 
is 

kT 1 d£ 

E c {r) = — jtt-t- 6 
3. Surface Model 

We use a simple, continuous model to include the effects of the surface on the conductivity. The 
fact that there are no carriers beyond the wire radius is accounted for by taking the volume 
density of carriers to decay spatially as a Gaussian with a width w beyond a certain radius r . 
A Gaussian function is chosen because of its mathematical simplicity and because its qualitative 
form is physically reasonable - we anticipate qualitatively valid physical conclusions using this form. 
This corresponds to a physical field that is proportional to the difference r — r for r > r and is 
directed towards the centre. An additional term proportional to 1/r in the effective internal field 
E c arises on using equation (5) due to the inclusion of the volume factor 2nr in the definition of /, 
as mentioned in section |21 

w( P )«,)H -:s exp( _ i) ;;:j m 




r > r ) 
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Figure 2. Radial profile of the relaxation time r(r). 



The relative roughness of the surface is described by lower relaxation times r(r) past a radius R. 
The change is also modelled by a Gaussian fall-off, with a width w r . For simplicity, the radii R and 
ro are assumed to coincide, and also the widths w and w r . The radial profile of the relaxation time 
is shown in Fig. El 

We take w to be of the order of 0.1 x R. 
4. Conductivity 

In the absence of a magnetic field, it is possible to make certain predictions analytically about the 
behaviour of the conductivity. We work in the limit where E is small, so that the drift velocity (v z ) 
is always small in comparison to the thermal mean velocity vth- If t{t) were a constant, tq, the 
drift velocity would be given by the familiar result 

(vz) = (10) 

m 

When t acquires a non-trivial radial profile r(r), its overall scale determines to what extent thin 
cylindrical shells at different radii affect their neighbours because of the diffusion of carriers from 
layer to layer. If this scale is so small that the mean free path of particles with velocities of 
the order of Vth is much smaller than the length scale ( over which r changes significantly (i.e. 
T i r ) v th C = w )i the spatial connectivity is negligible, and one may treat the different layers 
separately. In this case, one can define an effective relaxation time r geom which depends only on 
the geometrical distribution of the relaxation time weighted by the relative carrier concentration 
£(r). The use of £(r), which is the spatial factor in the unperturbed distribution function f , here 
is consistent with the fact that the relaxation-time-scale is small. 

fdr£(r)r(r) 
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For the model discussed in section |HJ this can be expressed simply in terms of the ratio u = w/R, 
providing a useful test for the numerical scheme in the limit of small r . 



It is therefore useful to consider the essential surface effect as the departure of the actual drift 
velocity for the nanowire as a whole from the value derived by inserting T geom in place of To in (3). 
This, as we have just observed, will be perceptible only when r v t h is comparable to the surface 
width and significant only when comparable to the wire radius. 

Now, it is important to note that the effect of spatially varying r in a conductor as a whole 
is asymmetric in its action, between smooth (high r) and rough (low r) regions. A rough region 
is much less affected by the gradient of r(r) in its neighbourhood than a smooth region. This is 
precisely because spatial connectivity is enhanced where r(r)v t h is large as explained previously. 
Thus, in a conductor which has both regions of high r and low r, the latter separately exhibit 
geometric (unconnnected) behaviour, showing no effect of the presence of the former. 

But the effect of spatial variation on regions with high r is to decrease mobility there, since 
carriers there can move enough to sample a significantly rougher region. Therefore, the effect of 
spatial variation on rough regions being negligible, its effect on a conductor as a whole, too, is a 
decrease in mobility. This consideration will later appear prominently in explaining the magneto- 
conductive effect as well (section EJ). Fig. El shows the variation of the conductivity as a function 
of the relaxation-time-scale for different temperatures. The ratio of the surface width to the wire 
radius u is fixed at 0.2. As r — > 0, the conductivity tends to the same fraction of the bulk value at 
all temperatures, which is seen to be almost exactly equal to the geometric factor 0.8815, the value 
obtained by putting u = 0.2 in equation (12). At higher temperatures, the departure from this ratio 
is also higher. This is physically expected since a higher temperature makes available higher radial 
velocities to the carriers in their random motion between collisions, thus increasing the communica- 
tion between different layers. Thus we see that even if the relaxation time were independent of the 
temperature, the conductivity would have a temperature-dependence in the low field limit because 
of the surface effect. Note that to isolate the surface effect, the complex temperature variation of 
relaxation time in a real medium is deliberately suppressed, though it can be included easily in the 
computational scheme. 

Further, it is seen in the inset that the data can be collapsed on to a single curve by using 



the transformation To — > tq^jT/Tq, where To is an arbitrary temperature. In other words, when 
B — the temperature is a reducible parameter in the limit of a low electric field, which can be 
accounted for exactly by rescaling r in proportion to the corresponding mean thermal velocity. 
This can be seen analytically by direct use of the Boltzmann equation, through a transformation to 
dimensionless variables. Rewriting eqn. (2) with the functional forms in (7), (8) and (9) we have 



1 + y/lTUJ + UJ 2 



(12) 



77 



geom 7~0 





+ qEr{r) 




p r d(j) qEr{r)p z 
m Or mkT 




(13) 
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Figure 3. Relative conductivity as a function of the relaxation-time scale To for different 
temperatures. er is the bulk value of the conductivity corresponding to a constant relaxation 
time To- The inset shows the data collapse resulting upon scaling to by a factor y/T/To. The fixed 
parameters are E — 26kV/m, R — 200nm and w = 40nm. 



Now introducing the variables u = p/VmkT and s = r/R, we get 

When E is small, <fi is of the same order of smallness, and therefore the second term in the left 
hand side of (14) is negligible. This leaves us with a linear differential operator containing r only 
in the combination t ^/T acting on <fi in the left hand side. The appearance of To and T in the right 
hand side (inhomogeneous term) of course, merely alters the overall scale of <fi. Hence the function 
0(u, r) may be written as -g=^ x H( u, s; tq\/T) where H is the solution of (14) with the constant 
-^|= absent in the right hand side. The drift velocity is 

' kT J d 3 u f ds u z (p(u, s) 



m j d 3 u J ds C,(s)e " 2 / 2 
Eqro J d 3 u J ds u z H(u, s; t ^/T) 
m J d 3 u j ds £(s)e~" 2 / 2 
/ d 3 u J ds u z H(u, s; ToVT) 



(v z )o 



fd 3 ufds£(s)e- u2 / 2 
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Figure 4. Relative conductivity as a function of wire radius. Here E — 26kV/m, T = 300K and 
w = 20nm. 

Thus (v z ) I (v z )o (= (x/cr ) depends on r only through the product ToVT. 

Fig. |3] shows the variation of the conductivity as a fraction of the corresponding bulk value with 
wire radius, parametrised by r . As r decreases, the results approach geometrically determinable 
values; further, as u increases, not only does the geometrical factor move further below 1 but also 
the actual conductivity departs more and more from the geometrical value. Expectably, if the width 
is fixed, increasing the wire radius diminishes the strength of the surface effect. 

The conformity with physical expectation and analytically known limits in the results above 
indicates the reliability of the numerical scheme and its suitability for other calculations based on 
the Boltzmann equation. 

5. Limit of Zero Surface Width 

A simple parametrization of the surface in a zero- width surface model is provided by the "specu- 
larity coefficient" e, which is the non-zero probability of carriers incident on the surface suffering 
a scattering event there: they undergo diffuse rather than specular reflection [46J. This is often 
encountered in the literature; for instance, it has also been used in the analogous context of thermal 
transport by phonons HS1 |4"9~1 15U|. We now consider the physically important limit of w — > 
which should relate the parameters in our model to the parameter e. Fig. El shows the numerical 
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Figure 5. Variation of conductivity with surface width for different relaxation-time scales. The 
geometric factor is included for reference. Here E = 26kV/m, T = 300K and R — 200nm. 

results obtained by varying the surface width parameter, along with the graph of the analytical 
geometric factor. Significantly, it is seen that even for large tq, the departure from the bulk value 
seems to go down to a small number, possibly zero. This is in contrast to pQ where the magnitude 
of the surface effect is unbounded as a function of the relaxation-time scale or mean free path. It 
may be noted that the computation becomes increasingly more expensive, as the surface width is 
decreased; we believe that the lowest width used here is low enough to allow us to draw qualitative 
conclusions from the results. 

We thus see an apparent contradiction between the results obtained with two different char- 
acterisations of the surface, which demands closer scrutiny. A direct contradiction results only if e 
in the abrupt surface model is considered a free parameter assignable arbitrary values between 
and 1; one could interpret the results we have seen as placing an upper limit on e not far above 
zero. But this, on the other hand, would imply that such a parametrization is not very meaningful. 
Indeed, a value of e = 1 is assumed in [I]. 

The contrast between the two approaches becomes understandable when one recognizes the 
correlation between e and the fraction of carriers resident within the surface layer. In a finite width 
model, this fraction is always finite; it approaches zero only in the zero width limit, when e —>■ 0. 
Equivalently, a non-zero value of e is only compatible with the existence of a finite fraction of car- 
riers on the surface. Thus if a finite width model is accepted as more realistic, a value of e close to 
1 is physically possible only when a large part of the total number of carriers resides on the surface 
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There is another way of pointing out the essential difference between an abrupt surface model 
and a finite width one. In the former, it is assumed that a scattering at the surface necessarily 
returns the carrier into the bulk of the surface; whereas this is not the case in the latter, once again 
bringing to mind the existence of a finite surface population. Finally, in the former, there is no 
scale of the relaxation-time at which the surface effect saturates, whereas, in the finite width model, 
saturation occurs where the mean free path exceeds a few times the length scale £. Saturation is 
physically expected in the latter case - if in some case, typical carriers can sample regions with 
significantly different values of r during a lifetime, any further increase of the lifetime should not 
make a big difference. 

A future investigation of the surface effect resulting from a carrier distribution involving non- 
trivial features near the surface rather than the simple, monotonic decay studied here is expected to 
provide further insight. In particular, it may provide an explanation for the apparent experimental 
evidence [HI] that e is close to 1. In any case, the considerations stated in this section challenge 
the conception of the abrupt surface model, which ought to be re-evaluated in the light of these 
numerical results. 

6. Longitudinal Magneto-conductance 

The effect of a longitudinal magnetic field (B) on the conductivity of a large nanowire is of par- 
ticular interest because the surface provides a classical kinetic mechanism for a magneto-resistive 
effect. In the large field limit, one would expect the magnetic field to decouple different spatial 
regions by confining the carriers kinematically, constraining them to helical motion between colli- 
sions. Thus, when the cyclotron radius is smaller than (, the carriers are effectively localised by 
the magnetic field, so that the conductivity is largely determined by the geometric factor. As we 
have seen in section HJ the effect of spatial connectivity in a wire with a rough surface is to lower 
the conductivity below the geometric factor. Since a large B field ought to reduce connectivity, the 
magneto-conductance clearly ought to be positive. Further, this magneto-conductance is expected 
to saturate with B when the cyclotron radius for carriers with typical momentum values (around 
the thermal momentum p t ^) becomes much smaller than (. On the other hand, the magnetic field 
has an appreciable effect only when the reciprocal of the cyclotron frequency for carriers with typical 
momenta is smaller than the relaxation time; the magnetic field is ineffective if the probability for 
a collision before the carrier velocity turns round once is high. 

At intermediate values of B, it is not obvious whether there can be a case of negative magneto- 
conductance; the field could conceivably tend to produce a net movement of carriers in some regions 
down the gradient of r(r). The results presented in Figs. |H1 and [7| suggest that this does not occur; 
the magneto-conductance is seen to be positive within the accuracy of the calculation in all cases. 
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Figure 6. Fractional longitudinal magneto-conductance as a function of relaxation-time scale for 
high magnetic fields. E = 26kV/m, T — 30QK, R = 200nm and w = 40nm. 

For relaxation times corresponding to the bulk mobility of a common semi-conducting material 
like Si, the effect of the magnetic field is surprisingly small. In particular, the saturation values of 
B exceed realistic laboratory values, whereas the magnitude of the fractional magneto-conductance 
is still only a few per cent. Again, this is in contrast to the results in pQ. Although it is to be noted 
that the latter pertain to the degenerate limit (metallic case), the crux of the difference lies in the 
surface model, as argued in sectional In short, if the surface effect is small, the magneto-conductive 
effect must also be limited by the corresponding departure from the geometric factor. 

However, there still remains a qualitative similarity between our results and the results in pQ; 
for instance, we see that the low B results lend themselves to a good parabolic fit, which is con- 
sistent with pp. Further, the arguments we have presented are substantiated by the fact that pQ 
systematically predicts a greater surface effect than the experimental results used for comparison, 
a discrepancy noted there itself. Since the discrepancy is despite the use of best-fit values, the 
parametrization itself must be regarded as dubious. It is important to note, however, that the 
results in pQ and our results are both consistent with recent experimental findings [23 E3 EI] to 
the extent that this surface magneto-resistive effect can be experimentally delineated from intrinsic 
quantum-mechanical effects. 

The intriguing fact that the magnetic field increases conduction in all cases examined prompts 
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Figure 7. Variation of fractional magneto-conductance with B going up to very large values. Fixed 
parameters used same as in Fig. [5] 



a second look at the general effect of a spatially varying relaxation time. For one thing, using the 
same numerical scheme, we find another curious fact - the magneto-conductance does not flip sign 
when the profile of r(r) is inverted, i.e. with a thin region of high mobility surrounding a core of 
low mobility, even in the high B limit. More generally, we could not find a profile of r(r) such as 
to yield a negative magneto-conductance. 

We attribute the absence of a negative magneto-conductance to the asymmetry between smooth 
(high r) and rough (low r) regions in the effect of a spatial variation of r, discussed in section HI 
As explained there, the surface effect, which is the result of connectivity in the model at hand, 
tends to decrease the conductivity below the geometric value. Thus the introduction of a magnetic 
field, whose basic action is to undo the transport connection between different regions by confining 
carriers, in reversing the effect of connectivity, can only increase the conductivity. 

One also sees, by the same line of reasoning, that the surface effect must be small if r(r) is 
assumed to be continuous, for, the region that contributes most to it is that where both r and its 
gradient are large; but this region is of the order of the surface width w, which is expected to be 
small compared to the wire radius. 
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7. Conclusion 

A direct solution of the Boltzmann transport equation offers a powerful approach to transport 
calculations in large nanowires. We have presented a simple finite width model in a cylindrical 
geometry and shown that despite its simplicity, its scope is significantly greater than that of an 
abrupt surface model, from which it exhibits qualitative differences. This challenges the utility of an 
abrupt surface model, especially when parametrized by the specularity coefficient. Our numerical 
results show that the classical magneto-conductive effect in large nanowires of materials like Si is 
limited to a few per cent even with magnetic fields beyond the range of a practical laboratory setup. 
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